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Abstract. - We find that three nonstoichiometric cubic binary chromium chalcogenides, namely 
Cr3S4, Cr3Se4, and Cr3Te4, are stable half-metallic ferromagnets with wide half- metallic gaps 
on the basis of systematic state-of-the-arts first-principles calculations. We optimize their struc- 
tures, and then calculate their magnetic moments, electronic structures, formation heats, and 
elastic moduli and investigate their structural stability and robustness of ferromagnetism against 
antiferromagnetic fluctuations. Our calculated results show that the three sulvanite phases are 
structurally stable and ferromagnetically robust, and hence could be realized as epitaxial thin 
films. We attribute the structural and ferromagnetic stability and the better half-metallicity to 
their special effective Cr valence 2.667+. These findings will open doors for much more high- 
performance spintronic materials compatible with current semiconductor technology. 
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■ It is believed that next-generation high-performance 
qomputers can be achieved through using the spin free- 
dom of electron in key materials and devices of current 
semiconductor technology [1,2]. The half-metallic fer- 
romagnet, first discovered in NiMnSb by de Groot et 
ql in 1983, has almost 100% spin polarization near the 
Fermi level [3-5] . This feature of half-metallic ferromag- 
iietism makes carriers have high spin-polarization near 
the Fermi level and avoid some spin-related scattering 
processes that should exist otherwise. These are essen- 
Ijial to practical spintronic applications [1, 2, 5]. Half- 
metallic ferromagnetism has been found in many materi- 
als, such as Heusler alloys [3,4,6], transition- metal oxides 
[7-9], and even graphcne nanoribbons under electric field 
[10]. The half-metallic ferromagnetic materials compatible 
with semiconductor technology is believed to be promis- 
ing candidates for achieving more powerful computers. 
For this purpose, many zincblende transition-metal pnic- 
tides and chalcogenides [11-18], regularly-doped semicon- 
ductors [19,20], semiconductor super-structures [21,22], 
and more related compounds have been proved to be 
half- metallic. Although great advance has been achieved, 
better half-metallic ferromagnets compatible with cur- 
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rent semiconductor technology arc highly desirable for the 
next-generation computers. 

Here we find three stable half-metallic ferromagnets, 
namely Cr3S4, Cr3Se4, and Cr3Te4, among nonstoichio- 
metric cubic binary transition-metal chalcogenides on the 
basis of state-of-the-arts first-principles calculations. We 
investigate their structural stability and ferromagnetic ro- 
bustness against antiferromagnetic fluctuations. Our cal- 
culated results show that the three nonstoichiometric sul- 
vanite phases have formation heats of down to -0.302 eV 
per formula unit with respect to corresponding zincblende 
phases and their half-metallic gaps can be up to 1.05 eV. 
We attribute the structural and ferromagnetic stability 
and the better half-metallicity to their special effective Cr 
chemical valence 2.667+, compared with 2+ in the stoi- 
chiometric structures. These nonstoichiometric cubic ma- 
terials have much better features for spintronic applica- 
tions based on semiconductors. More detailed results will 
be presented in the following. 

We use the full-potential linearized augmented-plane- 
waves method within the density functional theory (DFT) 
[23], as implemented in package WIEN2k [24,25]. The 
generalized gradient approximation (GGA) [26] is used for 
the exchange-correlation potential. Full relativistic effects 
are calculated with the Dirac equations for core states, 
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Fig. 1: (color online). The crystal structures of zincblende CrX 
(a) and sulvanite Cr3X4 (b) piiases, where X can be S, Se, or 
Te. The red (or black) ball denotes a Cr atom and the yellow 
(or gray) an X atom. 



and the scalar relativistic approximation is used for va- 
lence states. The spin-orbit coupling is neglected because 
it has little effect on our results. We use 2000 k points 
in the first Brillouin zone, make the harmonic expansion 
up to ?max— 10 in atomic spheres, and set i?mt x ^''^max to 
7.5. The radii of the atomic spheres of Cr and others are 
chosen so that as high accuracy as possible is obtained. 
The volumes are optimized in terms of total energy, and 
the internal position parameters with a force standard of 
3 mRy/a.u. The simplest antiferromagnetic structures are 
constructed by doubling the unit cells along the [100] and 
[110] directions. There are six Cr atoms in each of these 
doubled cells. We let three of the six Cr spins orient up 
and the other three down. All the spin values of the Cr, 
S, Se, and Te atoms are determined naturally by the self- 
consistent calculations. The structural stability is inves- 
tigated through deforming the structures (with the vol- 
umes fixed) along the three directions. The elastic moduli 
are calculated with the standard method implemented in 
WIEN2k [24,25]. The formation heats of the cubic sul- 
vanite structures are calculated with respect to the crys- 
talline Cr phase and the corresponding zincblende struc- 
tures. The self-consistent calculations are considered to 
be converged only when the absolute integrated charge- 
density difference per formula unit between the two suc- 
cessive loops is less than 0.0001|e|, where e is the electron 
charge. 

The ground-state phases for three stoichiometric CrX 
{X=S, Se, and Te) compounds are of nickel-arsenide (na) 
structure, but the zincblende (zb) structures of the CrX 
have been shown to be higher only by 0.28-0.36 eV per 
formula unit than the corresponding na ones [15], and 
therefore have very good stability. This actually stim- 
ulated experimental synthesis of zb-CrTe epitaxial thin 
films with a thickness up to 100 nm [16]. The sulvanite 
structure of Cr3X4 still has cubic symmetry with space 
group No. 215. Each of the cations has four anionic neigh- 
bors, but each of the anions three cationic neighbors. With 
the cations being at the (5, 0,0) sites, the anions occupy 
the {zx,zx,zx) sites. The two cubic structures are shown 
in Fig. 1. For all the sulvanite Cr3X4 compounds, the 



Table 1; The lattice constants a, the internal atomic parame- 
ters Zx, and the Cr-X bond lengths Z„ of the three sulvanite 
Cr3Jf4 phases. Presented in parentheses are those of the cor- 
responding zb-CrX phases. 



Name 


Cr3S4 (CrS) 


Cr3Se4 (CrSe) 


Cr3Tc4 (CrTc) 


a (A) 
Zx 

iv (A) 


5.344 (5.469) 
0.250 (0.25) 
2.315 (2.368) 


5.679 (5.833) 
0.250 (0.25) 
2.460 (2.526) 


6.117 (6.292) 
0.260 (0.25) 
2.686 (2.725) 



Fig. 2: (color online). Spin-dependent densities of states 
(DOSs) of the three sulvanite Cr^X^ phases for X = S (a), 
Se (b), and Te (c). The upper part in each of the three panels 
is the majority-spin DOS and the lower part the minority spin. 
The solid thick line represents the total DOS; the red (or gray) 
dash, blue (or light gray) dot, and black thin lines describe 
the partial DOSs projected in Cr and X atom spheres and the 
interstitial region, respectively. 




Stable half-metallic ferromagnetism in chromium chalcogenides 



lattice constant a and the internal structural parameter 
zx are optimized fully in terms of usual total energy and 
force standards, and then the Gr-X bond lengths 1^ are 
calculated with the optimized structures. We present our 
calculated a, zx^ and results for all the three sulvanitc 
Cr3X4 in Table [TJ Those of the corresponding zb CrX 
phases are shown in parentheses for comparison. Cr3S4 
and Cr3Se4 almost keep the same zx=0.25 as zb-CrS and 
zb-CrSe, and Cr3Te4 has zx=Q-260, a little larger than 
that of zb-CrTc. The sulvanite structures have smaller lat- 
tice constants than the zb ones by 1.6-2.6% and smaller 
Ci-X bond lengths by 1.4-2.6%. In the following, spin- 
dependent densities of states (DOSs), energy bands, and 
charge and moment density distributions are calculated in 
terms of the optimized structures. 
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(c) 

Fig. 3: Spin-dependent energy band structures (EBs) of the 
three sulvanite ^3X4 phases for X = S (a), Se (b), and Te 
(c). The left part in each of the three cases is the majority- 
spin bands and the right part the minority spin. The solid line 
with dots describes the EB structure along the high-symmetry 
points and the dot diameter is proportional to the Cr d weight 
at that point. 



The spin-dependent DOSs of the three sulvanite Gr3X4 
(X=S, Se, and Te) phases are presented in Fig. 2. For the 
convenience of description, wc label the four different sets 
of the energy bands as A, B, C, and D, respectively. For 
each of the three X cases, there is a narrow gap between 
A and B in the majority-spin (MAS) channel and a wide 
gap between C and D in the minority-spin (MIS). It is 
clear that there is a gap across the Fermi level in the MIS 
channel for each of all the three cases and therefore all 
the three arc half-metallic fcrromagncts. In addition to 
the total DOS, partial DOSs projected in the muffin-tin 
spheres of Cr and X atoms and in the interstitial region 
are presented too. It is clear that the B and D bands arc 
originated from Cr-d states and the A and C bands mainly 
have the X-p (S, Se, or Te) character. 

Presented in Fig. 3 are the spin-dependent energy 
bands, according to the corresponding DOS plots in Fig. 
2, for AT^S, Se, and Te. There are twelve bands in both 
A and C. They result from the fact that we have four X 
(S, Se, or Te) atoms in the sulvanite Cr3A4 unit cell and 
each of the four has three p orbitals. There are fifteen 
bands in both B and D because we have three Cr atoms 
here and each of the three has five d orbitals. The Cg 
bands are lower than the t2g ones. The Cr-s state is put 
to higher energy than the B and D bands in both of the 
MAS and MIS channels. This is in contrast with zb-CrA" 
phases whose Cr-s states are merged with the Cr-d t2g 
states. We have also calculated charge and spin densities 
in three typical planes. Our results show that the charge 
density from the A and C bands is distributed mainly in 
the neighborhoods of the X and Cr atoms, but the spin 
density is almost limited to Cr and X atoms only. The 
filled part of the B bands are substantially localized near 
Cr and X atoms. When X changes from S to Te, the X- 
p character becomes less and less in the electron density 
from the filled B bands because the ionic radius becomes 
larger and larger. 

Table 2: The magnetic moments (Af), half-metallic gaps 
(Ghm), magnetic energy differences (A_b), formation heats 
(-ffporm), and the elastic moduli {B, C' , and C44) of the three 



Cr3X4 phases, 


compared with the zb-CrX ph 


ises (in paren- 


theses). 








Name 


Cr3S4(CrS) 


Cr3Se4(CrSe) 


Cr3Te4(CrTe 


M (plb) 


10 (4) 


10 (4) 


10 (4) 


Ghm (cV) 


0.50 (0.07) 


0.83 (0.61) 


1.05 (1.00) 


Ab (eV) 


0.57 


0.56 


0.11 


-ffporm (cV) 


-0.148 


-0.114 


-0.302 


B (GPa) 


72.9 (63.2) 


57.7 (59.5) 


44.1 (45.9) 


C" (GPa) 


8.8 


5.9 (5.6) 


4.8 (5.5) 


C44 (GPa) 


34.2 


35.9 (50.7) 


32.9 (36.4) 



We summarize in Table [2] the total moments per for- 
mula unit (M), half- metallic gaps (Ghm) [27], and the 
energy differences per formula unit (A^; = Eaf — Efm) 
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for all the three sulvanite Cr^X^ phases (X=S, Se, and 
Te), where Eaf and EpM are the total energies for the 
AF and FM spin configurations. Here we define Ghm as 
the smaller of Ect and Eyt , where Ect is the bottom energy 
of the minority-spin conduction bands with respect to the 
Fermi level and Eyt the absolute value of the top energy of 
the minority-spin valence bands. Our GGA Ghm values 
should be rough estimates for the minimal energies for spin 
flip excitations. We obtain 10 fiB for the magnetic moment 
per formula unit in all the three cases. This is because we 
have ten Cr d electrons to contribute to the moment after 
eight of the eighteen Cr electrons of the three Cr atoms 
are bonding with the sixteen X p electrons of the four X 
atoms. The half-metallic gaps are larger than the corre- 
sponding zb-CrX phases because the Fermi levels move 
toward lower energy level due to the fact that there are 
one less Cr atom for CT3X4, compared with the cubic zb 
unit cell including four Cr and X atoms. The A^; results, 
from 0.11 to 0.57 eV, show that our FM solutions are very 
robust against possible AF fluctuations. 

In addition, we present the formation heats per formula 
unit (i?Form), the elastic moduH {B, C, and C44) of the 
the three Cr3X4 phases in Table [2] We calculate the for- 
mation heats to investigate the stability of the Cr3Ar4 with 
respect to the zb-CrX phases. For this reason, ifporm is 
defined as E{Ci3X4)-4E{zh-CTX) + E{Ci:), where E{f) 
is the total energy of the formula /. The negative values 
of i?Form for all the three cases mean that the sulvanite 
CT3X4 phases can be synthesized more easily than the 
corresponding zb phases, since zincblende CrTe has been 
synthesized in the form of epitaxial thin films [16]. The B 
and G44 values of the sulvanite phases are similar to those 
of the corresponding zincblende phases. The key tetrag- 
onal moduli, C, are 8.8, 5.9, and 4.8 GPa for the three 
cases. The latter two C values are comparable with those 
of the corresponding zincblende phases. Therefore, these 
calculated results show that the three sulvanite phases are 
mechanically stable. 

The nonstoichiometry plays key roles in the sulvanite 
Cr3X4. A Cr atom in the zb-CrX contributes a moment 
of ifiB, but three Cr atoms yield 10fj,B, not 12/x_b, in the 
sulvanite Cr3X4 because two more Cr electrons are needed 
to fill all the p orbitals of four X atoms. Effectively, the 
average valence of a Cr atom is equivalent to |-|- (higher 
than 2+) and the enhanced Cr valence makes the forma- 
tion heat negative and hence enhances the stability of the 
Cr3X4 with respect to the zb-CrX. In zb-CrX, Cr-dt2g 
states are hybridized substantially with X p ones, and Cr- 
deg are quite isolated. For the Cr3X4, Cr-d2;2_j,2 states 
remain isolated, and the other Cr d states are hybridized 
more strongly with the X p ones, which makes Cr s state 
higher in energy. 

In summary, we show that the three nonstoichiomet- 
ric cubic binary chromium chalcogenides are stable half- 
metallic ferromagnets with wide half-metallic gaps on 
the basis of systematic state-of-the-arts DFT calculations. 
Our calculated results indicate that the sulvanite Cr3S4, 



Cr3Se4, and Cr3Te4 are better than or approximately 
equal in both structural and ferromagnetic stability than 
zincblende CrS, CrSe, and CrTe, respectively. The struc- 
tural stability against crystal deformations and the nega- 
tive formation heats (-0.114~-0.302 cV per formula unit) 
mean that they can be synthesized, at least as epitaxial 
thin films, more easily than the corresponding zincblende 
phases. We attribute the structural and ferromagnetic 
stability and the better half-metallicity to the higher Cr 
chemical valence, 2.667-I-, in the sulvanite structures than 
2-|- in the corresponding zincblende ones. These findings 
will open doors for much more high-performance spin- 
tronic materials compatible with current semiconductor 
technology. 
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